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COMPUTER PREDICTION OF UHF BROADCAST SERVICE AREAS 

J.H. Causebrook, B.Sc. 



1. Introduction 

Current planning for UHF television requires recom- 
mendations for many relay stations to serve communities 
of more than 1000 people. Decisions on the location, 
power, aerial height etc. for these new transmitters may be 
determined by installing temporary equipment and measur- 
ing the fields so produced. This is a costly process, in both 
time and money, making it desirable to have a method for 
predicting the coverage of proposed transmitters by calcu- 
lation. Methods for making reasonably accurate predictions 
inevitably involve long and awkward computations and 
much data handling; thus manual calculations are not 
practicable where many cases have to be dealt with. Con- 
sequently computer programs have been written to perform 
these tasks. The purpose of this report is to describe the 
programs written to make field-strength predictions from 
small relay stations up to a maximum distance of 25 kilo- 
metres. For the longer distances needed for the calculation 
of co-channel interference, other programs have been 
written. 

It is not the intention in this report to give a full des- 
cription of all the operations performed in the computer 
program but simply to give an outline of the processes 
involved. 



2. The computer 

The programs are written for use on a terminal con- 
nected to a Univac 1108 computer. The programming 
language is Univac 1108 Fortran V but the differences 
from standard Fortran are not great. Therefore, with slight 
modifications, the program could be run on most com- 
puters. 



The first line of the input data to the field strength 
calculation program consists of the channel number, 
national grid reference of transmitter, receiving aerial height 
above ground level (if zero the aerial is assumed to be 10 
metres above ground level). An instruction 'noplot' may 
be written if only tabulated output is required (the form of 
the output will be discussed later). The second line of 
input data is any title required by the user and this is 
followed by several lines giving details of the transmitting 
aerial. 



4. The controlling sub-program 

The control ling sub-program is given the name 'CODER'; 
a flow diagram of its operations is given in Fig. 1. After 
the initial line of data has been read in, subroutines 'GR ID' 
and 'NGRGEO' determine the difference between true and 
grid north because aerial data are given relative to true 
north and profile data relative to grid north. The next 
significant operation is the use of a subroutine named 
'STADAT' which controls the reading of aerial data. The 
calculation of effective radiated power (e.r.p.) is discussed 
in Section 7. 

'CODER' now has the task of reading the topo- 
graphical data and is programmed so that if an error in 
format occurs the program stops. Thus the possibility of 
getting into an infinite loop is precluded. The subroutine 
'PROFTC turns the profile information into a series of 
heights spaced at regular intervals of 50 metres. Details of 
this process are given in Section 5. 

The subroutine 'CLUTTA' converts the information 
about buildings and trees, (i.e. where each group begins and 
ends) into a fractional coverage by the obstacles in each 50 
metre interval. 



3. The input data 

The topographical data are taken from maps and 
recorded on punched paper tapes using a machine called a 
'Terrain Coder'. An early version of this machine is des- 
cribed in BBC Research Department Report RA-22. In 
addition to the topographical data the Terrain Coder 
supplies information about the points at which field 
strengths are to be calculated. This may be either a set of 
individual points or a range of distances over which com- 
putations are made at 0-1 kilometre intervals. 

At present the Terrain Coder produces five-track 
paper tape in a code unsuitable for direct use by the 
program and a preliminary program named 'TERCOD' has 
been written to convert these data into an acceptable form. 
'TERCOD' also checks for certain errors which may have 
been introduced at the terrain coding stage. A listing of 
the output from 'TERCOD' may be obtained for manual 
checking and correcting before the information is used as 
input to the main program. 



Details of the calculations performed by the field 
strength subroutine 'TCFS' are covered in Section 6. This 
subroutine is used by 'CODER' once for each field point on 
a particular profile. When all the field calculations are 
complete a subroutine named 'OUTPUT is called; this is 
described in Section 8. 

If the first line of data has 'noplot' in the position 
allocated, the work of the program is complete. However, 
the user may require a plot, as shown in Fig. 2, giving the 
ffelds calculated along each of the profile lines from a trans- 
mitter, and this is achieved using subroutines contained in 
an available soft-ware package. The plot is made to scale 
such that it may be overlayed onto a map and the service 
area of the proposed transmitter quickly determined. 



5. Obtaining a good profile 

The data from the terrain coder contains contour 
heights and distances along a chosen profile line, together 
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with distances to the beginning and ends of stretches of 
water. At least two contours beyond the receiver point 
must be supplied so that interpolation can take place. 
These data are converted into heights at a regular series of 
distances spaced 50 metres apart. (This is done by a sub- 
routine named 'PROFTC which employs a further sub- 
routine named 'QUAD' to perform quadratic interpolation.) 
This process is essential at the end of the profile because of 
the way loss is calculated for buildings and trees and also 
provides heights between widely-spaced contours. The 
distances finally obtained are the points at which field 
calculations are made. 

A requirement, of the computer method for inter- 
polation, is that it should produce a profile as nearly equi- 
valent as possible to one drawn by hand. In addition the 
process should give a profile independent of its direction of 
production. The method adopted is described although it 
must be observed that the choice was partly subjective. 

At regular distances, D,-, from a terminal the heights, 
H- are derived from the weighting mean of two parabolic 
arcs. Let {h ^d^ ), [h^dj, (h^d^), and {h^d^) be four con- 
secutive contour height and distance pairs such that/): lies 
between d^ and d . A parabolic arc may be made to pass 
through ih^d^), [h^d^] and (h^d^) with constants a^, b^ 
and c^ and another through (h^d^), (h^d^) and (h^d^) 
with constants ag, bg and Cg. Thus at distance D-. two 
interpolated heights may be determined. 



(d) If the vertices of both parabolas lie inside the interval 
and are in the same sense, the interpolation is arranged 
to fall inside the triangle formed by the tangent to 
curve A at (h^d^), the tangent to curve B at (h^^d^ ), 
and the line joining the points [h^d^) and [h d ). 



The slopes of the tangents are: 



Sa = 2aAG?3 + b^ 



Sb = 2a^^3 + ^B 



(7) 
(8) 



The interpolated profile should produce a single peak 

and as in case (b) the slope is constructed without 

discontinuities. 

the junction of the tangent lines 



If D is less than the distance d^ of 



H-rK^ 



Sa + '/= 



,^3 - d^ 



-s 




(D. 



2 ' 

(9) 



and if D- is greater than d^ 



"r^^ 



Sg + yH 




(e) A limit is put onto the value of H- such that it shall 
not exceed the next contour va\ue to be crossed. 



^Bi = aB^)- +bgZ)j+Cg 



(1) 



The following rules are then applied to derive a 
weighted mean height//: from the above values//^: and 

^^Bi- 

(a) Water between d, and d , or d is the transmitter. 

\"/ 1 2 2 



^i=^Bi 

Water between d^ and d^ 



(2) 



(3) 



Water between d and d , or d is the transmitter. 

3 4' 3 



Hr^A] 



(4) 



(b) If the vertex of either parabola lies outside the interval 
d.d it is desirable for the profile to be made up of 
a curve with no discontinuity in slope. This is 
achieved with the interpolation: 



//; 



//AjW3-/^j)+//gj(/3j-6f^ 



(5) 



(c) If the vertices of both parabolas lie inside the interval 
d d and are in opposite sense the choice of the 
interpolation method is arbitrary and the simplest 
method was adopted i.e. 



6. The attenuation relative to free space field 

The principal subroutine for making calculations of 
attenuation relative to the free space field is named 'TCFS'. 
A flow chart of the procedures carried out by 'TCFS' is 
given in Fig. 3. 

The information entering the subroutine is mainly 
that describing the profile. The heights are not given as 
values above sea level but with reference to a cartesian 
frame formed by the vertical at the transmitter and the 
tangent to sea level at this point. In converting the heights 
an effective earth radius of four thirds the physical value is 
assumed. 

The first operation performed is to find the points 
which would be touched if a string were stretched between 
the transmitting and receiving aerials. A by-product of this 
process is the determination of the angle from the trans- 
mitter to its horizon or to the receiver whichwer is the less. 
This angle is needed to obtain the direction most relevant 
for calculating the e.r.p. 

The subroutine estimates the attenuation resulting 
from buildings and trees. This requires several applications 
of a function sub-program named 'FZ' which estimates the 
density of obstacles in a volume close to the receiver. The 
method of obtaining the attenuation due to buildings and 
trees is given in the Appendix. 



H: 



/^Ai+^Bi 



(6) 



The subsequent procedure depends on the answer to 
the question 'How many points, on the profile, are touched 
by the imaginary string?' Thus is the answer is: 
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Fig. 3 - Simplified flowchart for 'TCFS' 
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(a) 



'None' 



Subroutine 'DBMAXK' determines, from the profile, 
the points giving the maximum knife-edge diffraction loss. 
If this loss is insignificant the work of 'TCFS' is complete. 
However, if the loss is significant two forms of diffraction 
loss are employed: one the knife-edge loss and the other 
the diffraction loss due to a wedge-shaped obstacle. 

The method of calculating the loss over a wedge- 
shaped obstacle is desc-Jbed in BBC Research Department 
Report RA-91. In the current program however, the 

wedge surface is not defined as in Reference 3. The 
sty! isation of the profile is now made as follows: the apex 
is situated at the terrain point which gives the maximum 
knife-edge diffraction loss. The plane surfaces of the 
wedge are made to touch the terrain on either side of the 
apex, so as to make the internal wedge angle a maximum, 
provided this does not reduce the heights of the terminals 
above these surfaces to less than 15 feet. The part of the 
terrain used excludes that which is within the distances dj 
and d[^ kilometres from the apex, on the transmitter and 
receiver sides respectively. These distances are defined by: 



method described in Reference 3. To obtain a stylised 
surface in this case the first procedure adopted is to create 
a virtual edge from the terminal horizon rays, and then test 
the position of this edge in relation to the outermost edges. 
If the virtual edge and either of the outermost edges fall 
within the grouping criteria given by Equations (13) and 
(14), the virtual edge is used as the apex of the wedge 
sty I isation. 

If the virtual edge is not close to either of the outer- 
most edges, a surface is constructed which consists of four 
right-circular cylinders with slope continuity at their 'joins' 
(Reference 3). The outer cylindrical surfaces are made to 
fit the terrain in a similar way to that used for the wedge 
surfaces. 

In the multi-edge branch of the calculation, subroutine 
'DBMAXK' finds the edge which produces the maximum 
loss and the losses due to subsidiary knife-edges, as required. 
The approximation to the first term of the residue series, 
needed to calculate the diffraction loss over a cylindrical 
surface is performed with a function sub-program named 
'F'. 



dl = MAX jO-5, 



0-0994 {D- d^)d^ 
D + 0-0994 d. 



and 



d]^ = max] 0-5, 



0-0994 {D - d^]d^ 
1-0994Z)- 0-0994 d^. 



(11) 



(12) 



where/) = total path distance (km) 

(ig = distance from receiver to wedge apex (km) 

(b) 'One' 

The calculations are made for a single knife-edge and 
for wedge diffraction, with the apex at the point touched 
by the imaginary string. 

(c) 'More than one' 

The edges are examined to see if they may be put into 
groups. Groups are formed if: 



d 



i + 1 



d-< 0-5 km 



0-9095 < 



d;{D-d-,^,) 
d-^ 1 (D^d;) 



(13) 



(14) 



where d-^ and d-^ + , are a pair of adjacent edges touched by 
the string. Should it be necessary to group a set of edges, 
the set is replaced by a single virtual edge formed by using 
the first and last edges of the original set. This grouping 
may produce one virtual edge, in which case the calculation 
is made for a knife-edge and the wedge, but with d^_ and 
d^ spaced from the horizon points rather than the wedge 
apex. 

If two or more edges still exist a calculation is made 
for a series of knife-edges using the modified Deygout 



When the above operations are complete the program 
control is returned to 'CODER' with a series of values of 
attenuation relative to free space. These are Aqi_, A^^, A^ 
and A^ for clutter, knife-edges, cylinders and wedges 
respectively. 



7. The calculation of effective radiated power 

The e.r.p. from the transmitting aerial is required to 
determine the free space field at the receiver. The infor- 
mation needed for the calculation of free-space field is the 
transmitter power, details of the transmitting aerial, and the 
distance from transmitter to receiver. 

Two principal sub-programs process the aerial data. 
The first, 'STADAT', reads in information concerning the 
elements of which the aerial is composed, and the feeder 
loss. Information about the aerial may be given in the 
form of a tabulation of horizontal radiation pattern and 
vertical radiation pattern in which case a subroutine named 
'TABIN' is used. If the aerial pattern has the so-called 
'cardioid' shape it is not necessary to tabulate the pattern, 
for this is held in a subroutine named 'CARD16'. 

It is common practice for standard panels of 2X to be 
used in various combinations to form a transmitting aerial. 
In this case it is only necessary to state the orientation and 
relative position of each panel on the mast, the pattern 
being derived from a subroutine named 'UHPAN'. 

The program will also accept combinations of types 
of aerial elements. 

'STADAT' puts the h.r.p. and v.r.p. (both amplitudes 
and phases) for each element type into a 'common block' 
and returns the program control to 'CODER'. The flow 
chart, for the operations performed in 'STADAT', is shown 
in Fig. 4. 



(RA-125) 



COMMON BLOCK 
CONTAINING 
WAVELENGTH 



CONVERTS 
WAVELENGTH 
TO FREQUENCY 



CALL'CARD16' 



CALL'UHPAN' 




(RA-125) 



Fig. 4 - Simplified flow chart for 'STADA T 
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The second sub-program necessary to produce the 
e.r.p. is named 'AEPAT'. Its use is controlled from a sub- 
routine named 'OUTPUT' to be described in the next 
section. The main tasl< of 'AEPAT' is to compute the 
formula. 



k " 



i[Qk+Pk + 



27r 



— ( |X|^ -1- my 1^ -1- nZ|^)] 



(15) 



7. The combined attenuation relative to a free space 
field 

8. The effective radiated power in the direction of 
receiver 

9. The angle between the horizontal and the 'string' at 
the transmitter. The sign convention is positive 
downwards. 



where E is the field intensity in a direction defined by the 
direction cosines (I, m, n). Ttie relative drive current in 
the k* element is given by ^|^e' k and its radiation field 
in the specified direction is given by i?|^e' k. The effective 
position of the k* element in cartesian co-ordinates is 
given by (X|^,y^, z,^). 

The production of the various parameters required 
above is aided by two further subroutines named 'RELTP' 
and 'I NT'. Their use can be seen in Fig. 5 which is a flow 
ch art of the operations of 'AEPAT'. I n any one appi ication 
of 'AEPAT' 200 values of e.r.p. may be determined for 
various angles of elevation with one value of azimuth 
corresponding to a profile line. In practice the e.r.p. 
values produced are for a standard transmitter power of one 
kilowatt 



9. The optimisation process 

The values obtained in 'TCFS' are derived on a 
theoretical basis. However, to obtain such values simpli- 
fying stylisations of the profile are made and this means 
the figures cannot be relied upon without some form of 
empiricism. For this purpose a measurement project was 
undertaken. The measurements were within 50 metres of 
radial lines from the transmitters. The radial lines were 
divided up into 100 metre intervals and all measurements 
within the 100 metre square were used to produce a single 
value to represent the area. These were then taken as the 
values to be used in optimising the calculation system. The 
total number of such areas obtained was 377, these being 
along 17 radial lines from four transmitters. By this 
method it was hoped to cover as many types of path as was 
economically feasible. 



a The subroutine 'OUTPUT' 

As already stated the subroutine 'OUTPUT' controls the 
use of 'AEPAT' thus acquiring the e.r.p. values. The 
various losses produced by 'TCFS' are combined to produce 
a single attenuation value relative to free space; this com- 
bination process is described in Section 9. The field is 
produced using these values in the following form: 



F= 106-9- 20log^^D + P+ T-A 



(16) 



where F is the field strength in dB relative to 1 /uV/m (dB/u), 
D is the distance from transmitter to receiver in km, P the 
e.r.p. in decibels relative to one kW (for Tx power 1 kW), 
T is the transmitter power in decibels relative 1 kW and A is 
the attenuation relative to free space. 



The minimising of differences between calculation 
and measurement was done by a non-linear optimision 
technique based on a method by Dickinson. This involved 
programs which were used in the development stage only. 
They consisted of three subroutines; one to read in data 
and print results; another to compute the function to be 
minimised; and the third which did the search for a mini- 
mum. The final attenuation relative to free space was made 
up from the values given by 'TCFS' in the following 
manner: 



Xi^k + Xj^c + Xj^w + ^4 + '^s^CL + Xs^CL 



(17) 



where the various values of A are defined at the end of 
Section 6 and where x j etc. are values to be determined by 
optimisation. The function to be minimised was 



The subroutine 'OUTPUT prints results in tabular 
form. An example of a table for one profile line is shown 
in Fig. 6. The information at the head of the table is self- 
explanatory. The columns have the following meanings: 

1. Distance from transmitter to test receiver point. 

2. The predicted field-strength in dB/U, assuming the 
receiving aerial is at 10 metres above ground level, 
although other heights may be taken. 

3. The predicted loss as a result of trees and buildings on 
the path. 

4. The predicted loss assuming profile is a series of knife- 
edges. 

5. The predicted loss assuming profile is of cylindrical 
form. 

6. The predicted loss assuming profile is of wedge form. 



y = S(^-^» 



(18) 



where the summation is over all the available samples and 
Af^ is the measured attenuation relative to free space. 



The result obtained was: 



A = 0-503^1^ + 0-4bA^ + 0-3a4^ + 2-4 + 



+ 6-54 



CL" 



0-1754 



GL 



(19) 



The standard deviation of the differences between ^(yi and 
the values of A given by Equation (19) was 4-28 dB and the 
largest differences were 11-8 and —12-9 dB. These results 
are plotted in Fig. 7 and compared with the zero difference 
line. 
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Fig. 6- An example of the tabulated output 

10. Profile plotting on an x-y incremental plotter. Thus the user has a pictorial 

representation of his problem from which he may make 
A program has been written to plot the ground profile decisions such as altering certain transmitter parameters, 
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Fig. 7- Comparison of measured and calculated attenuations 



including more receiving sites whien situations are critical 
or removing them where results are obvious. He may 
decide that more profiles are necessary when he sees from 
existing profiles that the situation is critical in certain areas. 
The profile plot may also help to identify errors occurring 
in the terrain coding or interpolation processes. An 
example of a plot is given in Fig. 8 (the original is in two 
colours; red and black). The crosses along the profile line 
represent the original figures supplied by the terrain coder; 
the line is that produced by the interpolation using these 
values. Each receiving aerial position investigated is marked 
by a small cross on the top of a vertical line. Along the 
axis of the profile are representations of buildings and trees 
(trees being shown taller). These are not shown on the 
profile itself because this would present technical difficul- 
ties. 



11. Conclusions 

A computer program has been developed to predict 
service areas of relay stations as an aid to service planning 



in the UHF Bands. This program can be employed in 
suitable cases as an alternative to making measurements in 
the field which are relatively costly and time-consuming. 
As the average cost for obtaining the service area of a 
small relay station by means of the program is small, it 
will be appreciated that considerable savings may be made. 
In addition to the saving in cost, the prediction can be made 
in a relatively short time and therefore also produces a 
saving in effort. This program, in conjunction with other 
programs designed to take account of co-channel inter- 
ference, is already in use and experience is rapidly being 
gained in their operation. The accuracy obtained as com- 
pared with measurements appears to be adequate. 
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Appendix 
The Effects of Buildings and Trees 



Most domestic reception of broadcast signals takes 
place with receiving aerials at roof top level and for pre- 
diction purposes this is taken to be 10 metres a.g.l. At 
such a height, and at u.h.f., there is liable to be considerable 
attenuation caused by trees and buildings in front of the 
receiving aerials. Measurements have indicated that the 
average loss is about 10 dB but that it could be higher than 
30 dB. A routine has been devised to give a value for the 
loss due to this cause. 

An example of the problem requiring solution is de- 
picted in Fig. 9 which shows the receiver end of a path and 
demonstrates how the direct ray may be obstructed by 
buildings and trees. 

Ground heights are usually slowly varying (rarely 
does a slope exceed 45 degrees) but the buildings and trees 
have predominately vertical sides, are of many different 
heights and many of them usually occur in a small area. 
Because of the many paths that must be dealt with it is not 
possible for each obstruction to be identified and even if 
such data were collected it is doubtful if a method could be 
devised to use it. 

The most economic way of collecting the required 
data is from Ordnance Survey maps. However, this has 
several disadvantages: 

1. No heights of the objects are given. 



2. Only woodland trees are shown, not trees in hedge- 
rows, road-sides or gardens. 

3. The areas covered by buildings are exaggerated to a 
highly variable extent. 

The exaggeration of buildings shown on various types 
of maps was investigated and it was found that 1:63,360 
Ordnance Survey maps exaggerated by an average of 3 
times, with a large scatter about this value. The 1:25,000 
old series Ordnance Survey maps exaggerated by 1-5 times, 
but the new series is not significantly exaggerated. 

Without having definitive data on the heights of the 
required objects the next best approach is to employ some 
form of statistical information. A study of this subject 
resulted in a probability distribution for these objects 
exceeding a given height above ground. These are shown 
in Fig. 10 as fg and f-p for buildings and trees respectively. 

There is a high probability that the final signal will 
result from a complex combination of many forms of 
propagation mechanism. 

There are: 

1. Multiple diffraction over obstacles. 

2. Multiple diffraction round them. 

3. Many reflections off many objects (scattering). 

4. Transmission through a combination of obstacles. 




0-2 0-3 0-4 0-5 0-6 0-7 0- 

distance, km 
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Fig. 10 - Probability of obstacles exceeding a given height 

The theory for any of these is difficult; to combine 
them is not feasible. These means of propagation may 
produce a high location variation, because interference 
effects give the so called 'standing wave' patterns. Time 
variations may also result if for example these patterns are 
varied by the movement of trees in the wind. To reproduce 
these patterns by calculation is clearly not possible, but if 
some form of clutter data are given it should be possible to 
calculate a good representative figure for the area. 

Theories only exist for cases where features at right 
angles to the propagation path vary slowly or the changes 
are regular, as with a diffraction grating. However, it is 
possible to devise a method which takes some recognition 
of the various theories and at least gives an indication of 



the important parameters. It is then necessary to find a 
relationship between these parameters and a set of measure- 
ments. Ideal controlled experiments to determine para- 
meters by measurement can be thought of but it is not 
practical to adopt these. 

In the absence of the ideal approach all that can be 
done is to take a set of measurements covering as much of 
the problem as possible and devise a mathematical model 
which Includes various constants determined by minimising 
the error between measurements and calculations. 

Along these lines Kinase has written of a method 
using a single parameter symbolised by T. The strict 
interpretation of his parameter is the percentage of 'clutter' 
occupying the volume of the first Fresnel zone. However, 
Kinase allows the parameter to have a less rigid meaning 
for practical purposes; i.e. the percentage, of the area in 
question, occupied by 'clutter' exceeding 10 metres above 
ground level. Even with thisreduced meaning some explicit 
knowledge of heights is required. To obtain this would be 
too expensive for the purpose in hand. In any case F is not 
considered to be the best that can be done in terms of a 
parameter because the field received is highly dependent on 
the way in which the first Fresnel zone is filled. To illus- 
trate this point consider an absorbant screen placed in one 
of three positions (1, 2 and 3) shown in Fig. 11. The 
fractional volume filled will be the same in each case but 
the effect on the received field will be very different. 

Buildings will only transmit a small amount of energy 
so the principal mechanisms will be diffraction and reflec- 
tion. For simple knife-edge diffraction the relevant para- 
meter is the Fresnel 'v' variable. Little can be done about 
reflections except to say that when the diffraction by 
buildings severely reduces the direct signal, reflected fields 
will tend to increase the average received signal. 

Most practical cases will be covered if we consider a 
volume of revolution about the line of an assumed string 
stretched between transmitter and receiver over the ground 
profile with radius given by: 



\s{D- s) 
2D 



(20) 



where D is transmitter to receiver distance and x is the dis- 
tance from the receiver for the required radius. The v para- 
meter may be written in terms of this variable and the 
height, h, of the edge from the centre of the circle, i.e.: 



(21) 



However, in this particular problem there is no definitive 
knowledge of/;. The proposed procedures are based on the 
division of the volume (defined above) at right angles to the 
main axis is shown by the dotted lines in Fig. 11. Then 
values, F^, for the fractional filling of each of these slices 
may be derived using the formula: 
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Fig. 1 1 - Obstructions in the first Fresnel zone 



H + rf 



■nrt 



(H-y) 



iQ{y)dy (22) 



H 



where H is the height of the ray line above ground level, j 
is a variable height above ground and ^^ty) is the proba- 
bility of any object exceeding heightj. 

Within the context of the problem in hand, trees 
differ from buildings in two respects; firstly they are 
usually taller and are thus more liable to be above the top 
of the volume and secondly they allow a signal to be trans- 
mitted through them which may well exceed the signal 
diffracted round or over them. A paper on the subject of 
this transmitted energy was published by Saxton and 
Lane where results are given in terms of decibels per 
metre. If this were the only mechanism acting, the loss (in 
decibels) along a ray line would be given by: 



D 



p(s) 6s 



(23) 







where p{s) represents the likely density distribution of 
trees over a reasonable distance and x is the loss per unit 
distance through the trees. Of course, x may be a function 
of several factors such as frequency, polarisation, type of 
trees and season. 

The mechanism of Equation (23) will predominate 
when the diffracted signal is low and the depth of trees is 
not great as is likely to be the case with trees close to the 
receiver. Under these conditions the loss through a length 
AS may be written as: 



^^S 



xF^AS 



(24) 



If many intervals have clutter in them an attempt 
must be made to combine these, but there does not 



appear to be any feasible analytic way of solving this 
problem. However, knife-edge diffraction theory does lead 
to two conclusions which can be used. Firstly, the overall 
diffraction loss is dominated by the greatest of the indi- 
vidual losses and secondly this domination is the greater 
the closer subsidiary edges are to the dominant edge. For 
buildings this means it is reasonable to take only a maxi- 
mum value for each small interval, and then emphasise the 
maximum value of many of these small intervals by taking 
the square root of the sum of squares to these values. 

The work of Saxton and Lane and Sofaer and Bell ' 
leads to a figure of 0-2 dB per metre at u.h.f., for attenua- 
tion through trees, and ignoring differences due to polarisa- 
tion and season should not lead to errors of more than 
4 dB. However, for a certain depth of trees the largest 
signal will be received by diffraction over their tops. Thus 
if it assumed that the diffraction loss caused by a clump of 
trees at, say, 1 km from a receiver is about 10 dB diffraction 
will be the dominant mechanism if the trees are more than 
50 metres deep. If AS is chosen as 50 metres and each 
interval is multiplied by the thickness of trees in it the 
result will be in close agreement to the values of Reference 
6. With these thoughts in mind the loss parameter is taken 
as: 



iCL- 



2(Pbi^bi + Pti^ti)' 



(25) 



where i labels each interval of width AS 

Pbi is if no buildings are shown on the map in the 

interval and 1 if any buildings are shown 
F^j andF^i are assumed fillings of the volumes in 

question 
P^l is the fraction ground coverage by trees, shown 

on the map in the appropriate interval. 

In practice the main contribution to the parameter 
will come from the part of the path which is close to the 
receiver. This is because in the rest of the path the top of 
the volume is well clear of clutter, transmitters being built 
with aerials higher than those of domestic receivers. Errors 
will not usually be greater than a few decibels if the only 
clutter values used are those within 2 kilometres of the 
receiver. 
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The ground is likely to be included in the volume, ever, it is not the intention to calculate the field at a point 

but the loss resulting from this is dealt with separately so it but to obtain a representative field for a small area, thus 

is not reasonable to include it again. some 'averaging out' of errors is bound to occur. Nonethe- 

less, large errors may occur if a tall building lies in 

It is evident that the influence of buildings and trees front of the aerial, or many trees exist there and are not 

on -field strengths cannot be calculated definitively. How- shown on the maps. 
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